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METHOD AND SYSTEM FOR TIMING RECOVERY AND 
DELAY SPREAD ESTIMATION IN A COMMUNICATION SYSTEM 



TECHNICAL FIELD OF THE INVENTION 
The present invention relates generally to communication systems, and 
in particular, to a new timing recovery and delay-spread estimation scheme 
for communication systems that employ cyclically extended symbols. 

BACKGROUND OF THE INVENTION 
Timing recovery in communication systems is related to the process of 
identifying symbol boundaries in a received signal, so that each symbol can 
be windowed and processed separately so that its value can be determined. 
In some communication systems such as OFDM (other examples include 
IFDMA, Single Carrier with cyclic extension), a guard interval is inserted 
between successive symbols to overcome inter-symbol interference (ISI) 
caused by multipath delay-spread in the communication channel. Usually 
each symbol is cyclically extended with a prefix and/or a postfix to cover the 
guard interval. The cyclic extension absorbs the delay-spread and thus keeps 
the data portion of the symbol free of ISI. When the channel delay-spread is 
less than the duration of the cyclic extension, only a portion of the cyclic 
extension is corrupted while the rest remains ISI-free. This creates ambiguity 
in the timing recovery process because there is more than one possible 
position of the sampling-window Tor obtaining an ISI-free representation of the 
symbol. 

FIG. 1 shows an exemplary OFDM symbol with cyclic prefix partly 
corrupted by ISI and the range of ISI-free sampling-window positions available 
for that symbol. 

Although known methods for timing recovery in OFDM systems are 
adequate and beneficial in many situations, they present several 
shortcomings. In a method proposed by Jan-Jaap van de Beek et a/., "ML 
Estimation of Time and Frequency Offset In OFDM Systems", IEEE 
Transactions on Signal Processing, Vol. 45, No. 7, July 1997 (hereinafter 
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"Jan-Jaap van de Beek "), the cyclic extension is used to identify the best 
sampling position. Unfortunately, this method is optimized and well suited 
only for a single-ray channel. As a consequence, in a multipath-fading 
channel, the sampling position picked by this method may not be free of ISI 
5 and it can waver depending on which of the multipath rays is the strongest at 
a particular time. These shortcomings in a multipath-fading channel are 
critical, considering that one of the purposes behind OFDM is to reduce the 
effects of multipath fading. 

In another method proposed by T.M. Schmidl et a/., "Low-Overhead, 

10 Low-Complexity [Burst] Synchronization for OFDM", Proceedings of ICC 

1996, vol. 3, pp. 1301-1306 (hereinafter "Schmidl"), a special training symbol 
is used to estimate a sampling position. The disadvantage of this method is 
that the sampling position picked by this method could jump about within a set 
of valid positions, leading to jitter in the timing estimates. The jitter makes it 

15 difficult to use averaging, e.g., a phase locked loop (PLL), to obtain a steady 
sampling position. Another disadvantage with this method is that the timing 
estimate is based on the training symbol alone and there is no averaging over 
the fading process. This creates problems in a fast fading scenario where the 
instantaneous delay profile at the training symbol could be significantly 

20 different from the true delay profile. 

Both of the above-mentioned methods for timing recovery in OFDM 
systems involve finding the correlation between two symmetric portions of the 
OFDM symbol. In the Jan-Jaap van de Beek method, the symmetry between 
the cyclic extension and the data portion of the OFDM symbol is exploited, 

25 whereas in the Schmidl method, the symmetry between the first and second 
halves of a special training symbol is exploited. 

In addition, the Jan-Jaap van de Beek and Schmidl methods both 
employ a sliding sum that adds up consecutive correlation values over the 
window of symmetry. For example, in the Jan-Jaap van de Beek method, 

30 consecutive correlation values are summed together over a window of length 
L, where L is the length of the cyclic extension. Consequently, the correlation 
is summed over the length of the entire cyclic extension, regardless of the fact 
that part of it could be corrupted by ISI. This degrades the correlator output 
by adding incoherent correlation values to the coherent ones. It can also 
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result in a loss of information regarding how much of the cyclic extension is 
corrupted, or in other words, the delay-spread. 

Accordingly, there is a need for an improved method of timing recovery 
and delay-spread estimation in communication systems employing cyclic 
5 extension. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 illustrates an exemplary OFDM symbol; 

FIG. 2 is a flowchart illustrating a method of performing timing recovery 
10 in accordance with the invention; 

FIG. 3 is a block diagram of an exemplary receiver in accordance with 
the invention; 

FIG. 4 is a set of graphs showing an exemplary comparison between 
the ensemble correlation function of the present invention and the correlation 
15 functions of two known methods of timing recovery for OFDM; 

FIG. 5 is a set of histograms showing an exemplary comparison 
between the timing recovery method of present invention and two known 
methods of timing recovery for OFDM; 

FIG. 6 is a flowchart of a method for estimating channel delay-spread in 
20 accordance with the present invention; and 

FIG. 7 is a flowchart of a method of processing an information signal 
based on an ensemble correlation metric in accordance with the invention. 
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DETAILED DESCRIPTION OF THE 
PRESENTLY PREFERRED EMBODIMENT(S) 

The present invention provides a timing recovery scheme that 
demarcates a complete range of ISI-free sampling positions available in an 
OFDM symbol. The invention also features a method for estimating the 
multipath delay-spread in the channel. The delay-spread estimate provides a 
receiver valuable information about the nature of the channel. The receiver 
can use this information, for example, to adjust the frequency-domain channel 
interpolation filter bandwidth in order to improve performance for the given 
channel condition. 

The invention uses the correlation between the cyclic extension and 
the data part of the OFDM symbol to estimate the timing. However, it has key 
differences from known techniques in the way the correlation is computed and 
used to estimate timing. For example, according to an aspect of the 
invention, a correlation function provides better delay-spread and timing 
information by not combining together consecutive correlation values over the 
entire length of the cyclic extension. Instead, the function separately 
maintains the correlation value for each sample position within the OFDM 
symbol duration and combines the correlation values for identical sample 
positions over an ensemble of OFDM symbols. 

In this manner, the invention can provide an improved method for 
estimating the symbol timing from the received OFDM signal. Each sample 
position in the OFDM symbol can be inspected separately. Consequently, the 
method can provide a steady, jitter-free estimate of the symbol timing for a 
variety of channel conditions. In addition, the method can average the 
correlation function over a large number of symbols to provide robustness to 
fast fading, i.e., to limit the impact of instantaneous fades on performance. 

The correlation function retains delay-spread and timing information by 
not combining consecutive correlation values together over the entire length 
of the cyclic extension. Instead, it separately maintains the correlation value 
for each sample position within the OFDM symbol duration and combines the 
correlation values for identical sample positions over an ensemble of OFDM 
symbols. Identical sample positions in successive OFDM symbols are 





15 



20 



separated by A/+Z. samples, N being the fast Fourier transform (FFT) size and 
L being the length of the cyclic extension, where N and L are integers. 
Therefore, the correlation function sums correlation values spaced by A/+/_ 
positions. This sum is done over an ensemble of, for example, M OFDM 
symbols and then normalized, where M is an integer. If r denotes the 
received signal sequence, this "ensemble correlation function" can be 
expressed in equation form as follows: 



Since the ensemble correlation function computes the correlation 
coefficient separately for each sample position in the OFDM symbol duration, 
it gives information about what portion of the cyclic extension is corrupted due 
to channel delay-spread. This makes it possible to estimate the delay-spread 
and also to identify and demarcate a range of ISI-free sampling points 
available. In a single-ray channel, where there is no delay-spread, the 
ensemble correlation function has a plateau of width L, indicating that every 
sample of the cyclic extension matches the corresponding sample of the data 
portion of the symbol. In a multipath channel, the width of the plateau 
reduces depending on the extent of delay-spread D in the channel. 

FIG. 2 is a flowchart 20 of an exemplary method for using the 
ensemble correlation function of Eq. (1) to estimate OFDM symbol timing. In 
step 22, the ensemble correlation function is computed over a plurality of 
OFDM symbols according to Eq. (1 ). Next, in step 24, the output of the 
ensemble correlation function is post-processed as follows. First, a 3-tap or 
5-tap median filter can be used to remove any sharp glitches in the correlation 
function, while still retaining the distinct nature of the plateau and its sharp 



jV (* + m(N + L))r(k + m{N + L) + N) 




0<k<N+L 
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Next, in step 26, the peak value, p of the correlation function is found 
and a threshold value, t is determined as a function of the peak as t = p-a(1- 
pj, where a > 0 is a design parameter. 

In step 28, the first and last points where the correlation function 
crosses the threshold are determined and the region in between these two 
points is declared the ISI-free sampling region. The width of this region W is 
subtracted from the total length of the cyclic extension to obtain an estimate of 
the channel delay-spread. 

D = L-W (2) 



In the initial acquisition stage, the plateau may be split between the two 
ends of the ensemble correlation function. In other words, the plateau starts 
at a point close to k = A/+L, and wraps back around to k = 0. One way to 
alleviate this wrap-around is to circularly shift the correlation function such that 

15 the peak value p is positioned at the center, i.e., at k = (N+L)/2. This ensures 
that the plateau lies approximately at the center of the correlation function. An 
equivalent adjustment can also be made to the timing estimate as well. 

In a tracking mode, the ensemble correlation function is computed and 
processed as above for every group of M OFDM symbols received. The 

20 estimated start and end points of the plateau are tracked using two phase 
locked loops (PLLs). This yields a steady estimate of the delay spread and 
the ISI-free sampling region. 

The ensemble function of Eq. (1) is one example of a correlation 
function usable with the present invention. Other forms of ensemble functions 

25 can be used, for example, those having different normalization functions or 
denominator computations than those disclosed herein. 

FIG. 3 is a block diagram of an exemplary receiver 50 in accordance 
with another embodiment of the invention. The receiver 50 can be included in 
a communication system, such as a cellular and/or paging system, using an 

30 OFDM scheme or the like. The receiver 50 includes a timing recovery/delay- 
spread estimator 52, that includes a correlator 54, a filter 56, a peak (max) 
detector 58, a comparator 60, a pair of PLLs 62-64, and a subtractor 66. 
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The correlator 54 receives digital samples representing OFDM symbols 
and computes the ensemble correlation function in accordance with Eq. (1). 
The output of the correlator 54 is filtered by the filter 56, which can be any 
suitable digital filter, including the either of the median filters discussed above 
5 in connection with FIG. 2. The peak detector 58 determines the peak value 
and provides as output to the comparator the threshold value described above 
in connection with FIG. 2. 

The comparator 60 compares the filtered outputs of the ensemble 
correlation function to the threshold value to determine crossing points. The 
10 PLLs 62-64 can be any suitable PLL, digital or analog, for tracking the 
threshold crossings. 

The subtracter 66 outputs the delay-spread estimation by subtracting 
the width of the sampling region from the length of the cyclic extension of the 
OFDM symbol. 

15 The architecture of the timing recovery/delay-spread estimator 

(TRDSE) 52 is shown as an example. The structure of the TRDSE 52 can 
include any suitable combination of hardware and/or software components for 
performing the functions of the components 54-66 described herein. 

FIG. 4 shows a plot 84 of the ensemble correlation function used by the 

20 present invention and compares it to the correlation functions used by the 

known methods presented by Jan-Jaap van de Beek 80 and Schmidl 82. Two 
different channel conditions are shown in the columns of the plot - one where 
the channel is ideal with no delay-spread (left column) and another where the 
channel is noisy and has some delay-spread (right column). The number of 

25 OFDM symbols used in computing the ensemble correlation function shown 
was M = 200. It can be observed from the FIG. 4 that the method of Jan-Jaap 
van de Beek provides a single peak and does not give any delay-spread 
information. This peak does not always occur at an ISI-free sampling 
position. The method of Schmidl provides a plateau, but this plateau is not 

30 "well-marked" and is therefore difficult to identify and demarcate. 

In contrast, the method of the present invention provides a distinct 
plateau, which clearly demarcates a range of ISI-free sampling positions 
available. Thus, the timing information extracted by the present invention is 
superior compared to the previous methods. 
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FIG. 4 also illustrates that the total height of the plateau is proportional 
to the total signal power and the height of each individual step that leads to 
the plateau is proportional to the power of a corresponding multipath ray. 

FIG. 5 shows a histogram 104 of the timing estimate generated by the 
5 method of the present invention as compared with the known methods of 
Jan-Jaap van de Beek 100 and Schmidl 102. The parameters used to 
generate the histogram were N = 1024, L = 256, M = 200, a = 0.35. The 
OFDM symbol rate was 20 kHz. The channel was a 4-ray channel with a 
delay span equal to one-half the length of the cyclic extension. An 
10 exponentially decaying power delay profile was used with the last ray having a 
power 3 dB below the first ray. Each ray was faded at a Doppler frequency of 
300 Hz, which corresponds to a vehicle speed of 58 mph at a carrier 
frequency of 3.5 GHz. The signal-to-noise ratio (SNR) was 10 dB. 

It can be observed from the histogram plots 100-102 that the Jan-Jaap 
15 van de Beek method can frequently pick sampling positions outside the valid 
sampling region. Thus, the timing estimate from this method is not free of ISI. 

The Schmidl method picks sampling positions within the valid sampling 
region most of the time, but these estimates jump around within the valid 
sampling region leading to jitter in the timing estimate. 
20 In contrast, the method of the present invention provides an accurate 

estimate of the start and end of the valid sampling region. Any particular point 
within this region can then be fixed as the estimated sampling position. This 
method is much more reliable than previous methods in the sense that the 
timing estimate is ISI-free for a variety of channel conditions. 
25 Furthermore, since this method provides an estimate of the entire valid 

sampling region, it is possible to sample the signal using more than one 
sampling window and combine the outputs appropriately to improve 
performance. Also, the delay-spread information provided by this method can 
be used to adjust the frequency-domain channel interpolation filter bandwidth 
30 to provide optimum performance. These techniques can yield up to 1 dB 
improvement in performance. 
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While specific embodiments of the present invention have been shown 
and described, it will apparent to those skilled in the art that the disclosed 
invention may be modified in numerous ways and may assume many 
5 embodiments other than those specifically set out and described above. 
Accordingly, the scope of the invention is indicated in the appended claims, 
and all changes that come within the meaning and range of equivalents are 
intended to be embraced therein. 
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CLAIMS 

1 . A method for timing recovery in an orthogonal frequency division 
5 multiplexed (OFDM) system, comprising: 

computing an ensemble correlation function over a plurality of 
OFDM symbols; and 

comparing the output of the ensemble correlation function to a 
threshold to define a sampling region within the duration of an OFDM symbol. 

10 

2. The method of claim 1 , wherein the sampling region defines the 
boundaries of the OFDM symbol. 

3. The method of claim 1 , further comprising: 

15 filtering the output of the ensemble correlation function. 

4. The method of claim 3, further comprising: 

using a median filter to filter the output of the ensemble 
correlation functions. 

20 

5. The method of claim 1 , further comprising: 

determining a peak value included in the output of the ensemble 
correlation function. 

25 6. The method of claim 5, further comprising: 

determining the threshold as a function of the peak value. 

7. The method of claim 1 , wherein the threshold crossing points of 
the ensemble correlating function output define the sampling region. 

30 
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8. A system, comprising: 

a correlator for computing an ensemble correlation function over 
a plurality of received OFDM symbols; and 
5 a comparator, operatively coupled to the correlator, for 

comparing the output of the ensemble correlation function to a threshold to 
define a sampling region within the duration of an OFDM symbol. 

9. The system of claim 8, further comprising: 

10 a filter, operatively coupled to the correlator, for filtering the 

output of the ensemble correlation function. 

1 0. The system of claim 8, further comprising: 

a max detector, operatively coupled to the filter, for determining 
15 a peak value included in the output of the ensemble correlation function. 

1 1 . The system of claim 8, for use in a receiver. 

12. The system of claim 1 1 , wherein the receiver is wireless. 

20 

13. The system of claim 8, further comprising: 

at least one phase locked loop for tracking the edges of the 
sampling region. 



25 



11 



14. A method for estimating delay spread in an orthogonal 
frequency division multiplexed (OFDM) system, comprising: 

computing an ensemble correlation function over a plurality of 
5 OFDM symbols, each of the OFDM symbols including a cyclic extension 
having a length; 

comparing the output of the ensemble correlation function to a 
threshold to define a sampling region within the duration of an OFDM symbol; 
and 

10 subtracting a width of the sampling region from the length of the 

cyclic extension of the OFDM symbol to obtain an estimate of the delay 
spread. 

1 5. The method of claim 1 4, further comprising: 

15 filtering the output of the ensemble correlation function. 

1 6. The method of claim 1 5, further comprising: 

using a median filter to filter the output of the ensemble 
correlation functions. 

20 

17. The method of claim 14, further comprising: 

determining a peak value included in the output of the ensemble 
correlation function. 

25 1 8. The method of claim 1 7, further comprising: 

determining the threshold as a function of the peak value. 

19. The method of claim 14, wherein the threshold crossing points 
of the ensemble correlating function output define the sampling region. 

30 
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20. A system for estimating delay spread in an orthogonal frequency 
division multiplexed (OFDM) system, comprising: 

a correlator for computing an ensemble correlation function over 
a plurality of OFDM symbols, each of the OFDM symbols including a cyclic 
extension having a length; 

a comparator, operatively coupled to the correlator, for 
comparing the output of the ensemble correlation -function to a threshold to 
define a sampling region within the duration of an OFDM symbol; and 

a subtracter, operatively coupled to the comparator, for 
subtracting a width of the sampling region from the length of the cyclic 
extension of the OFDM symbol to obtain an estimate of the delay spread. 

21 . The system of claim 20, further comprising: 

a filter, operatively coupled to the correlator, for filtering the 
output of the ensemble correlation function. 

22. The system of claim 20, further comprising: 

a max detector, operatively coupled to the filter, for determining 
a peak value included in the output of the ensemble correlation function. 

23. The system of claim 20, for use in a receiver. 

24. The system of claim 23, wherein the receiver is wireless. 
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METHOD AND SYSTEM FOR TIMING RECOVERY AND DELAY 
SPREAD ESTIMATION IN A COMMUNICATION SYSTEM 

5 

ABSTRACT OF THE DISCLOSURE 
A timing recovery scheme demarcates a complete range of 
inter-symbol interference free (ISI-free) sampling points available in an 

10 orthogonal frequency division multiplex (OFDM) symbol. The timing recovery 
scheme computes a correlation between the cyclic extension and the data 
part of OFDM symbols to estimate timing. The correlation function retains 
delay-spread and timing information by not combining consecutive correlation 
values together. Instead, the function separately maintains the correlation 

15 value for each sample position within the OFDM symbol duration and 
combines the correlation values for identical sample positions over an 
ensemble of OFDM symbols. In this manner, the scheme not only provides 
timing estimates, but can also provide estimates of the multipath delay-spread 
in a channel. The delay-spread estimates provide valuable information about 

20 the nature of the channel. A receiver can use this information, for example, to 
adjust the frequency-domain channel interpolation filter bandwidth in order to 
improve performance for the given channel condition. 
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Compute a "correlation value" for each received 
sample by taking the conjugate product of the sample 
with another sample that is N samples away 
A(k) = r'(k)r(k + N) 



Compute an "energy value" for each received 
sample by computing the square of its magnitude 
B(k) = \r{kf 



J For each of the A/+L sample positions within one OFDM j 
| symbol period, compute the numerator of an "ensemble i 
j correlation metric" by combining the correlation values • 
! across an ensemble of M OFDM symbols 



^A(k + m(N + L)) 



For each of the N+L sample positions within one OFDM 
symbol period, compute the denominator of the ensemble 
correlation metric by combining the energy values across 
the ensemble of M OFDM symbols 



£B(k + N + m(N + L)) 



Compute the ensemble correlation metric, 
p(k) by dividing N(k) by D(k) for k=0 to N+L-1 

D(k) 



Median filter the ensemble correlation j 
metric to remove glitches 



Find the peak, p of the ensemble correlation 
metric and compute a threshold, t as a function 
of the peak value 



Find the plateau region of the ensemble 
correlation metric by determining the portion 
that exceeds the threshold. Declare this region 
of width, Wthe ISI-free sampling region. 



Compute an estimate of the channel delay-spread. £ 
by subtracting the width of the plateau, Wfrom the 
length of the cyclic extension, L 
D = L-W 
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FIG. 6 
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j For each of a plurality of sample positions within a time 
j span of a symbol block, compute an "ensemble 
I correlation metric", where the metric is dependent on 
I pairs of received sample values, and where a pair of 
! received sample values is comprised of the two 
I received sample values at times tO + NTs and to + NTs • 
j + T1 , respectively, where N is an integer, Ts is the 
j extended symbol block length, and T1 is the symbol 
I block length. 



£ 

Determine a sampling position for demodulating 
the received signal based on the plurality of 
ensemble correlation metrics. 
OR 

Determine a channel channel impulse response 
duration based on the plurality of ensemble 
correlation metrics. 
OR 

Modify the processing used to detect the 
information signal based on the plurality of 
ensemble correlation metrics. 
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